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ABSTRACT 

Because of the complex nature of radiative heat-transfer calculations 
in many practical cases , especially where surfaces are involved, the use of 
digital computing machines becomes mandatory. It follows that methods that 
take advantage of the assets of these computers should be utilized. The 
Monte Carlo method is discussed as a means of solving radiative-transfer 
problems, and the literature dealing with its application is reviewed. An 
example problem is presented and discussed with regard to a Monte Carlo 
solution. Applications of the method to general problems in radiation are 
pointed out, and the areas where advantages over more conventional techniques 
exist are examined. 
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CALCULATION OF RADIANT HEAT EXCHANGE 


BY THE MONTE CARLO METHOD 

By John R. Howell 

Lewis Research Center 
Cleveland, Ohio 

SUMMARY 

Because of the complex nature of radiative heat -transfer calculations 
in many practical cases, especially where surfaces are involved, the use of 
digital computing machines "becomes mandatory. It follows that methods that 
take advantage of the assets of these computers should he utilized. The 
Monte Carlo method is discussed as a means of solving radiative -transfer 
problems, and the literature dealing with its application is reviewed. An 
example problem is presented and discussed with regard to a Monte Carlo solu- 
tion. Applications of the method to general problems in radiation are pointed 
out, and the areas where advantages over more conventional techniques exist are 
examined. 

INTRODUCTION 

The transfer of energy by thermal radiation for all but extreme situations 
is considered to have a solid theoretical foundation. The chief difficulty in 
solving practical problems involving radiative transfer does not lie, therefore, 
in formulating the problem. Rather, the resulting formulation, especially when 
other modes of heat transfer must be considered, gives equations that are so 
complex in many cases that complete solution by analytical techniques is im- 
possible, and solution by standard numerical methods impractical. 
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This difficulty arises because of the nonlinearity of the energy balance 
equations. Here combined energy transfer is occurring by radiation, where 
fourth-power temperatures are present, and by conduction and/or convection, 
involving first -power temperatures. -Further, the radiative terms in these 
equations take the form of multiple integrals, while conduction terms contain 
second derivatives. Further still, the radiative surface properties which 
appear are often functions of wavelength, direction, and temperature. If 
gases are involved, these variables plus pressure can strongly affect the 
local gas -radiation properties. The complete energy balance on each element 
of the system then takes the form of a nonlinear integro-differential equation. 

Because of this, three approaches to radiant -interchange calculations have 
been used. First, radiation is assumed negligible, thus admitting defeat. 
Secondly, the problem is simplified to the point where solution is possible by 
the means at hand. This is done by making as many assumptions, reasonable if 
possible, as are necessary, and philosophically accepting the resulting loss 
of accuracy if not validity. Surfaces that are black, gray, diffuse, or spec- 
ular and gases that are opaque, transparent, gray, or isothermal are assumptions 
that fall into this category. Few problems in radiative transfer are solved 
without explicitly or implicitly making one or more of these assumptions. 
Finally, as a third choice, more advanced techniques can be developed in the 
hope of accurately handling the mathematics of the more complex situations. 

One step along this third path is the Monte Carlo technique. Monte Carlo 
is in essence a method for statistical sampling of events. This sampling is 
done by analysis of a model of the real physical problem being considered. The 
method was developed to a high degree by those working on the first atomic bomb. 
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In solving problems in thermal radiation, the method has some very great 
advantages over conventional methods. Of course, there are also a few dis- 
advantages . 

The Monte Carlo method, its application to thermal radiation problems, 
and its advantages and disadvantages for these problems will now be discussed 
at some length. 

THE MONTE CARLO METHOD 

Monte Carlo was developed as a way of treating problems in which the 
happenings at a given location are known, at least in the form of statistical 
distributions, but in which the equations that describe the interactions 
between locations are extremely difficult to solve. One example is the local 
neutron flux produced by the diffusion of neutrons in the core, reflector, 
and shielding of a nuclear reactor. Because the neutrons arising from the 
fission process undergo different sequences of absorption, fission, and 
scattering, a solution for the neutron flux at all points in a heterogeneous 
system can be difficult. However, the frequency of events occurring along 
the path of an individual neutron are fairly well understood. This leads to 
the idea of following sample neutrons, and determining the events along their 
paths by picking events at random from the appropriately weighted set of 
possibilities at each point. By letting each neutron sample represent a group 
of real neutrons, and by following enough samples, the flux at each point can 
be determined. 

Of course, such a large number of simple calculations must be performed 
that a digital computer becomes a necessity. However, the fact that Monte 
Carlo depends on a large number of simple repetitive calculations and decisions 
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means that it is ideally suited for the computer. 

APPLICATION TO THERMAL RADIATION 

Applying the Monte Carlo method to thermal radiative transfer problems 
involves setting up a physical model which characterizes radiative processes. 

In neutron diffusion, sample neutrons obeying known scattering, absorption, 
and fission laws meet this need. In radiation, the corollary sample particle 
is the photon. However, if the photon itself is chosen as a particle to 
sample, the problem arises that the wavelength of the photon depends on its 
energy. It is more convenient to choose a common energy for all samples to 
be followed. Using this technique, each sample then becomes a bundle of 
photons with the same wavelength. Each sample can have a different number 
of photons according to its wavelength. However, the total energy carried 
by every bundle is the same. 

With these bundles of energy as our samples, it becomes relatively 
straightforward to simulate radiative processes. 

For an example, examine the energy transfer between element dA-^ at 
temperature T^ and surface Ag, an infinite plane, at temperature T 2 = 0 
(see Fig. l). So that some previous statements about directional and spectral 
properties will gain substance, let element dA]_ have emissivity 

(l) 

and let area 2 have emissivity 

g 2 = ^ ) ( 2 ) 

and assume only that the emissivity of both surfaces is independent of 
circumferential angle G (Fig. l) . This is the case for real surfaces prepared 
by sandblasting, plating, or etching. 
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For such a surface, the total emitted energy per unit time is 


Q e,l " e 1,l aT l ^1 


( 3 ) 


where ^ is the total hemispherical emissivity given in this case by 

jt/2 

€ 1^*A 1 sin P cos P cLA 


e T = 



( 4 ) 


oT*i 


and i-^ is the Planck spectral distribution of blackbody radiant intensity. 
If it is assumed that Q , the total energy emitted per unit time by 

e > 

dA-L, is composed of N of the energy bundles described above , then the energy 
of each bundle , c, is simply 


c 


%,i 

N 


( 5 ) 


To determine the energy transferred from element dA-^ to surface A^, we now 
follow N bundles of energy through their emission from dA-^, and determine 
the number S£ absorbed at A^. If the energy reflected from A^ back to dA-^_ 
is neglected, the energy transferred per unit time from dA^ to A-p will be 

^ 1-2 “ cS 2 

The next question is how is each individual bundle path determined and 
how is a wavelength assigned to each bundle? However this is done, the directions 
and wavelengths of the N bundles must conform to the constraints given by the 
emissivity of the surface and the laws governing radiative processes. For 
example, if we assign wavelengths to N bundles, the spectral distribution of 
emitted energy generated by the Monte Carlo process (comprised of the cNA dA 
for discrete intervals dA) must closely approximate the spectrum of the actual 
emitted energy (plotted as jte-^i-^ dA versus A). To assure this, a number 
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of methods are available for choosing the energy-bundle properties (Refs. [1 
and 2]). Leave the radiation problem momentarily for a physical interpretation 
of two common methods of choosing sample properties. 

Choosing individual events . - Of the two methods to be outlined for 
randomly selecting events in a manner that obeys the physical constraints of 
the problem, the first is probably the most intuitively satisfying. It con- 
sists of choosing events directly from the curve of known probability of an 
event . 

Consider the probability distribution of Fig. 2(a), given by the relation 

.2 

pU) = — (7a) 

1000 v ' 


in the interval 0 < £ < 10 and p(|) =0 elsewhere. Normalizing this 
relation by the area under the curve of Fig. 2(a), gives 


P(s) 


= p(l) = 
"^ 1 Q p(S)« " 1000 


(7b) 


Such a normalized probability curve is called a probability density function. 

To choose values of | in such a manner that Eq. (7b) is satisfied by 
the distribution of chosen values, proceed as follows: 

Two numbers, R a and R^, are chosen at random from a large set of numbers 
evenly distributed in the range (0-1). On the digital computer, such R values 
are selected randomly by means of either a random number generating subroutine or 
from stored sequences of random numbers. Choosing R values at random ensures 
that each event in a history will be independent of preceding events. 

"tPhe two random numbers are then used to select a point ( p (|),|) on Fig. 2(b) 


by setting 
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P(|) = Ra; i = R b (| max - i min ) = 10R b 

This point is then compared to the value of P(£) at £ computed from Eq. (7b). 
If the randomly selected value lies above the computed value of P(£), then 
the randomly selected value of £ is rejected and two new random numbers are 
selected. Otherwise, the value £ that has been found is used. Referring again 
to Fig. 2(b), it is seen that such a procedure assures that the correct fraction 
of | values selected for use will lie in each increment d| after enough 
selections are made. 

The difficulty with such an event choosing procedure is that in some cases 
a large portion of the values of £ may be rejected. A more efficient method of 
choosing £ is therefore desirable. The method to be outlined is more efficient 

for many of the distributions which occur in radiative -transfer Monte Carlo 

calculations. 

This method is, in short, to integrate the probability density function 
P(|) using the relation 

R ~ P(l')d|' (8) 

^-00 

where R can only take on values in the range (0-1) because of the properties 
of P(£)* Eq. (8) is known as the cumulative distribution function. The 
function R is then taken to be a random number, and values of | are obtained 

by choosing R at random and solving Eq. (8) for the corresponding value of £. 

To show that the probability density of £ chosen in this way corresponds to 
the required P(£), we can again examine the probability density function of 
Fig. 2(b) . 
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Inserting the example P(|) of Eq. (7b) into Eq. (8) gives 

‘S 3 

1000 



R = / P(|’)d|' = 


O) 


Equation (9) is shown plotted in Fig. 3. Divide the range of f- into a 
number of equal increments A|. Suppose M values of R are now chosen in 
the range 0-1, and these M values are picked at equal intervals along R. 
There will be M values of | which correspond to these M values of R. 

The fraction of the M values of £ which occurs per given increment A| is 
then 



But AR/A| is of course an approximation to (dR/d|) if a large enough value 
is chosen for M and small increments A| are examined. But dR/d| can be 
seen from Eq. (7) to be simply P(|) and it has been shown, therefore, that by 
choosing values of f in this manner the required probability distribution is 
indeed satisfied. 

A similar procedure for use when the probability distribution is a non- 
separable function of more than one variable (for example, a reflectivity, 
which for real surfaces may have the angles of reflection (p, ©) depending 
on one another) is demonstrated in Ref. [1]. 

Selection of events for the example problem . - To return to the problem 
at hand, the wavelength of emission for the emitted bundle must be chosen. 

It is assumed here that the surface properties are product functions of 
the two variables angle and wavelength, that is 

e(A,p) = e(A)e(p) 
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This assumption is valid for many actual surfaces, since the wavelength vari- 
ation of emissivity, for example, rarely depends on angle of emission. It 
follows, therefore, that dependence on either variable may he found by integrating 
out the other variable. Then the normalized probability of emission occurring in 


the interval dA is 

»*/2 

p (A) = I F(A,|3)dp = 


2 ^r /2<5i <*' p)i * sin 3 cos 3 d3 dA 




(ii) 


Substituting into Eq. (7) gives 

l A/V 2 


2 it 


R A = 



e-^C A ' , 3)i^ i si n 3 cos 3 d3 dA' 


V*! 


( 12 ) 


If the number of bundles N is very large, and this equation were solved for A 
each time an R^ were chosen, computing time could become too large for practical 
calculations. To circumvent this problem equations like Eq. (12) can be numerically 
integrated once over the range of A, and a curve can be fitted to the result. 

A polynomial curve is often adequate, as is the case in this problem, giving 


A = A + BR a + CRjJ + 


(13) 


This equation rather than Eq. (12) is used in the problem-solving program. 

Following a similar procedure for the cone angle of emission 3 gives the 

relation /*3 Z' 

: J J e i(A,3')i^ sin 3' cos 3' dA d3 1 



2n 


P(3',A)dA d3' = 




which can again be curve -fit to give 


(14) 


3 = A' + B'Rp + C’R R + ... 


(15) 
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For a gray, diffuse surface, Eq. (12) reduces to 


R A,gray ~ 




1 V *A’ 


crT; 


= F 0-A 


(16) 


where Fq_-^ is the well known fraction of hlackhody emission in the wavelength 
interval (0 - A). Equation (14) for this case reduces to 

T 1 

R p,gray = 2 / sin p' cos p' dp' = sin 2 p_ L (17) 

Jo 


or, 

sin p-j_ = -y/R^ (18) 

The point to be made here is that computational difficulty in obtaining A 
from either Eq. (13) or (16) is not greatly different, nor is it much different 
for obtaining p^_ from either Eq. (15) or (18). The difference is mainly in the 
auxiliary numerical integrations of Eqs. (12) and (14), which are performed once to 
get the curve fits for the nongray-nondiffuse solution. As far as the problem- 
solving program is concerned, the more difficult case may just as well be solved. 
Thus increasing problem complexity leads to only gradual increases in the com- 
plexity of the Monte Carlo program and similar gradual increases in computer time. 

For emission of an individual energy bundle from surface dA^, then, a wave- 
length A can be chosen from Eq. (13), and a cone angle of emission Pp can be 
chosen from Eq. (15). There remains only specification of the circumferential 
angle 0p. Because of the assumption made earlier that emission did not depend 
on 0p, it is easily shown by the formalism outlined, and is also fairly obvious 
from intuition, that 0^ can be determined by 

6 1 = 2*R 0 (19) 
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Because the position of plane A 2 with respect to dA-^ is known it is a 
simple matter to determine whether a given energy bundle will strike Ag after 
leaving dA^ in direction (3q, 0q) • (It will hit Ag if cos 0q > 0, as shown 
in Pig. 1. ) If it misses , another bundle must be emitted from dAq. If the 
bundle strikes A^, it must be determined whether it is absorbed or reflected. 

To do this, geometry is used to find the angle of incidence p 2 of the bundle 
onto Ag. 

cos 32 = sin 3]_cos (20) 

Knowing the absorptivity of A 2 from Kirchoff 's Law 

= (^l) 

and having determined the wavelength A of the incident bundle from Eq. (13) 
and the incident angle 3 2 from Eq. (15), the probability of absorption of the 
bundle at A 2 can be determined. The probability of absorption is simply the 
absorptivity of A 2 evaluated at P 2 and A because the directional spectral 
absorptivity ag(A,P) is the fraction of energy incident on A 2 (in a given 
wavelength interval) from a given solid angle that is absorbed by the surface. 

The absorptivity is therefore the probability-density function for the 
absorption of incident energy. It is now easy to determine whether a given incident 
energy bundle is absorbed by using the first of the two event -choosing methods 
just outlined; that is, by comparing the surface absorptivity A 2 (A,p)> which 
corresponds to P(g) (the probability of absorption), with a random number 1^. 

If 

Ra< a 2U,0) (22) 

the bundle of energy is absorbed, and a counter Sr, is increased by one. Other- 
wise, the bundle is assumed to be reflected, and is henceforth neglected. 


This 
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neglect is reasonable if the absorptivity of Ag is large, or if the direc- 
tional reflectivity is such that few bundles are reflected back in their 
original direction. If not, angles of reflection must be chosen from known 
directional reflectivities, and the bundle followed further along its path 
until absorption or loss from the system. For the purposes of this example, 
little is to be gained by following the bundle after reflection from surface 
Ag because the derivation of the necessary relations is similar to that 
already presented. The bundles are therefore neglected. 

A new bundle is now chosen and its history followed. This procedure is 
continued until all N bundles have been emitted from dA 1 . The energy 
absorbed at Ag is then calculated from Eq. (6), 

This completes derivation of the equations needed for solution of the 
example problem. In putting together a flow chart (Fig. 4) , some shortcuts can 
be used. For example, the angle 0 can be computed first. If the bundle is 
not going to strike Ag on the basis of the calculated 0, there is no point 
in computing A and 3 for that bundle. Alternately, because 0 values 
are isotropically distributed, it can be noted that exactly half the bundles 
must strike Ag. Therefore the calculated 0 values can be constrained to 
the range < 0 < rt/2. If N bundles are emitted in this range, then 

the calculated heat transfer will be 

*1-2 = ^ (23) 

The solution of this problem by Monte Carlo is now complete. An 
astute observer will note that this example could be solved without much 
difficulty by standard methods. A more astute observer might note further 
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that extension to only slightly more difficult problems would cause serious 
consequences for the standard treatments (for example, consider introducing a 
third surface with directional properties into the problem and accounting for 
all interactions). On the other hand, the author has found few radiation 
problems that will not yield to ;a Monte Carlo, approach. 

Some problems that have been solved by Monte Carlo are now discussed. 

ADVANTAGES AND APPLICATIONS OF MONTE CARLO 

Several means are available in the literature for solving radiative transport 
problems. These methods have been called "standard" or "conventional" methods 
herein, and include the techniques developed by Poljak [3], Hottel [4, 5], 
Oppenheim [6], and Gebhart [7] as well as formulation in terms of integral 
equations. Each of these has advantages for certain types of problems, and 
all will outshine the Monte Carlo approach in speed and accuracy over some limited 
range of radiation calculations that is outlined roughly by the complexity of 
the problem. 

The chief usefulness of Monte Carlo to the thermal radiation analyst lies 
in this fact: Monte Carlo program complexity increases roughly in proportion 

to problem complexity for radiative interchange problems. This is an important 
advantage because conventional methods increase roughly with the square of 
complexity of the problem, due to the matrix form into which they fall. How- 
ever, because Monte Carlo is somewhat more difficult to apply to the simplest 
problems, it is most effective in problems where complex geometries and variable 
properties must be considered. In complex geometries, Monte Carlo has the 
advantage that simple relations will specify the path of a given energy bundle, 
whereas most other methods involve explicit or implicit integrations over 
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surface areas. Such integrations become difficult when a variety of skewed 
or curved surfaces are present. One application of Monte Carlo to such a 
problem was in Reference [ 8 ], where the directional reflectivity and absorp- 
tivity of conical cavities with diffuse internal reflections were calculated. 

Here a straightforward Monte Carlo technique gave an easily programmed solu- 
tion to a problem of not inconsiderable difficulty for conventional approaches. 

Monte Carlo was first applied in thermal radiation to problems involving 
gases. This had been suggested in Ref. [9], and the idea was applied in a 
series of papers (Refs. [10 to 12]) which treat radiation through gases in 
various geometries. One paper by Howell and Perlmutter Ref. [ 12] gives solu- 
tions for the temperature distribution and heat transfer through a nongray gas 
with temperature dependent properties, contained between infinite parallel black 
plates. The cases of plates with much different temperatures, and of plates with 
the same temperature and with a parabolic distribution of energy sources in the 
gas are solved. 

Fleck [13] about the same time presented a lengthy discussion of combined 
conduction and radiation in absorbing-emitting media using a Monte Carlo approach, 
including transient effects. 

Howell, Strite, and Renkel [14, 15] used Monte Carlo in the more practical 
problem of predicting the combined radiative and convective energy transfer in 
the rocket nozzle of a gaseous-core nuclear propulsion system. They allow 
^in Ref. [l5])the absorption coefficient of the propellant to vary with any two 
variables (from wavelength, static pressure, and static temperature) evaluated 
locally, and account for axial and (in Ref. [16]) radial temperature gradients. 

In addition, the ability of a layer of optically thick gas injected along the 
nozzle wall to attenuate the extreme predicted radiative fluxes was demonstrated 


in Ref . [ 16] . 
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Such a system, involving combined heat-transfer modes in a flowing 

system with variable properties and an arbitrary axisymmetric bounding 
surface would at best be a laborious test of the more conventional approaches 
to radiative transfer. Yet the Monte Carlo program used was able to predict 
temperature profiles in the propellant, converged to within 0.1 percent at 
points, along with axial heat-flux distributions at the nozzle wall. All 
this was done in less than 10 minutes of machine time. 

DISADVANTAGES OF MONTE CARLO TECHNIQUE 

Monte Carlo calculations give results which fluctuate around the real 
answer because the method is a repetitive experiment on a mathematical model 
used in place of the actual physical situation. The uncertainty can usually 
be found by applying the standard statistical tests, and can be reduced in the 
same manner as experimental error; that is, by averaging over more tests 
(bundle histories). 

No rigorous test exists to guarantee the convergence of Monte Carlo 
results to valid solutions. This has not as yet proven to be a difficulty in 
thermal radiation problems . It would often be immediately obvious that con- 
vergence to invalid solutions was occurring because of the limiting solutions 
and physical constraints which are known for most radiation problems. 

An example of the scatter of results and convergence criteria for a 
radiation problem is shown in Fig. 5, taken from Ref. [12] - This figure 
shows typical scatter in Monte Carlo calculations of the temperature dis- 
tribution in a real gas between diffuse black infinite parallel plates at 
different temperatures. The absorption coefficient k of the gas is a 
function of wavelength and local gas temperature. Because the limiting 
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diffusion and transparent solutions are known in this case and bracket the 
Monte Co.rlo results, the likelihood that convergence to an invalid solution 
is occurring is small. Also, the solutions follow what is known to he 
physically necessary; for example, the gas temperature increases monotonically 
from the colder to the hotter wall, and the gradient of emissive power in the 
gas varies in the expected manner with the property variations. 

The results in Pig. 5 are compared in Ref. [12] to a number of approximate 
solutions involving various simplifying assumptions for the gas absorption 
coefficient, and agreement within the limits of the assumptions is good. The 
results presented in Fig. 5 required an iterative procedure because of the 
dependence of the gas absorption coefficient on temperature. Convergence to 
within 2 percent required three iterations, each involving 10,000 bundle 
histories . 

Sample size ■ - Other difficulties can and do arise in sampling techniques, 
most of which have to do with obtaining an optimum sample size. For example, 
in the gas radiation problem of Fig. 5, the free path of an energy bundle 
between absorptions becomes quite short in optically thick gases. Any energy 
bundle undergoes a large number of absorptions and emissions in the gas before 
traversing the distance between parallel plates. This means that the computer 
time needed to obtain a sample of bundles transferred between the plates become 
prohibitively large. Conversely, for nearly transparent gases, relatively few 
absorptions occur in any given element of gas, and an insignificant energy 
sample is available to calculate the temperature of the gas element with 
accuracy. This accounts for the relatively greater scatter for the 0.5 centi- 

meter plate- spacing results of Fig. 5. 
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Both these problems are common enough in transport processes that are 
mathematically related to radiative transport that special methods of 
"weighting" the free paths of bundles have been developed (e.g., see Ref. [2]) 
to obtain adequate samples, thereby saving computer time and gaining accuracy, 
albeit at the expense of added complexity. 

In the gas-radiation problems available so far in the literature, it has 
been found that problems which cannot be accurately handled by either the trans- 
transparent or diffusion approximations fall into the range of mean free paths 
that allow efficient Monte Carlo programming without reverting to weighting 
techniques . 

CLOSING REMARKS 

Monte Carlo is discussed in the preceding section as a method suitable 
for use in the solution of complex problems in radiative transfer.. A sample 
problem is outlined to demonstrate its application, and some of the advantages 
and disadvantages of the technique are discussed along with pertinent literature 
references . 

From all this, certain conclusions emerge. First, Monte Carlo appears to 
have a definite advantage over other radiative -transfer calculation techniques 
when the difficulty of the problem being treated lies above some undefined level. 
Just where this level is cannot be established, probably being a function of the 
experience, competence, and prejudice of the individual working the problem. 
However, problems above this nebulous benchmark in complexity can be treated with 
greater flexibility, simplicity, and speed. Monte Carlo does lack a kind of gen- 
erality common to other approaches in that each problem may require an individual 
technique, and a dash of ingenuity often helps. This places a greater burden on 
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the programmer's backlog of experience and intuition where standard methods 
may allow programming through "cookbook" application of their formalism. 

Second, for thermal-radiation problems, the parameters and mathematical 
relations involved lie in ranges which allow straight-forward Monte Carlo 
programming without the need of the more exotic schemes occasionally necessary 
in other Monte Carlo transport studies. 

Third, with all its advantages, the method suffers from certain problems. 

The worst of these are the statistical nature of the results, and the lack of 
guaranteed convergence. It should be noted that the latter fault is common to 
all methods when complex problems are treated. 

Finally, it must be co mm ented that the person using Monte Carlo techniques 
often develops a physical grasp of the problems encountered, because the model 
being analyzed is simple, and the mathematics describing it are therefore on 
an unsophisticated basis easily related to the physical model This is in con- 
trast to the rather poor physical interpretations and predictions which we can 
make when working with, say, a matrix of nonlinear second-order integro- 
differential equations. 

LIST OF SYMBOLS 

A surface area 

c energy per sample bundle 

F_ , fraction of total blackbody emissive power that lies in the wavelength 

0-A 

range 0-A 

i, Planck spectral distribution of blackbody radiant intensity 

A 

(energy)/ (unit projected surface area) (unit solid angle) (unit wavelength 
interval) (unit time) 

p(i) 


probability distribution 
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P(| ) probability density function, Eq. (8b) 

Q energy per unit time 

R number chosen at random from set of numbers evenly distributed in range 

(0 - 1 ) 

S number of energy bundles absorbed at a surface per unit time 

T ab solute temperature 

a absorptivity of a surface 

P cone angle of emission (Fig. l) 

or cone angle of incidence 
e emissivity of a surface 

e T total hemispherical emissivity, Eq. (4) 

K absorption coefficient 

A wavelength 

i variable 

o Stefan-Boltzmann radiation constant 

0 circumferential angle of emission (Fig. l) 

Subscripts 
e emitted 

T total 

1,2 refers to surface dAp or A 2 , respectively 
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Fig. 4. - Computer flow chart for example program. 
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Fig. 5. -Emissive power distribution in hydrogen 
between parallel plates; no gas heat source; 
plate temperatures; T 0 = 9500° K; Tj = 4500° K 
(from Ref. 12 ). 











